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INTRODUCTION 

Many aircraft carry heavy concentrated weights attached 
to the wings, in the form of engines, auxiliary fuel tanks or 
weapons. These heavy weights can produce deliterious effects 
on the flutter speed of the wing. If they are not located at 
a propitious position, an undesirable lowering of the flutter 
speed results. Recently attempts have been made to adapt the 
concepts of active control to increase the flutter speed, but 
this method requires a complex system of electronics, sensors, 
controls, and activators. A passive method of achieving the 
desired effect would be more desirable from several standpoints, 
namely, less complex, less weight and more dependable. A 
passive system proposed by Mr. Wilmer Reed of NASA, Langley 
Research Center involves the concept of mounting the store on 
a pitch pivot having a very low pitch stiffness relative to 
the wing stiffness. It is the purpose of this research grant 
to analytically study the proposed concept. 

The first task was to investigate the concept utilizing 
a two-dimensional approach involving 4 degrees of freedom, 
namely, wing bending, wing torsion, store pitch and store 
vertical translation. This preliminary analysis was very 
encouraging and the results demonstrated that if the uncoupled 
store pitch frequency was below the wing bending frequency 
that the flutter speed was greatly increased. A second more 
complete analysis was developed utilizing a three dimensional 
structure, but retaining the two-dimensional, incompressible 
unsteady airforces of Theodor sen. The details of the analysis 
are given in the Appendix. 
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SYMBOLS 


a 

location of wing elastic axes, measured from 
wing midchord, based on wing half chord, 
positive aft 

b 

wing half chord 

El 

wing bending stiffness 

EMF 

-’"p/'^m 

EMS 


F,G 

real and imaginary parts of Theodorsen Function C(k) 

g(i) ,g{2) , 
g(3) ,g(4) 

structural damping coefficients for first four 
uncoupled bending modes 

g(5) ,g(6) , 
g(7) 

structural damping coefficients for first three 
uncoupled torsional modes 

g(8) ,g(9) , 
g(io) 

structural damping coefficients for store vertical 
translation, pitch, and fore and aft, respectively 

GJ 

wing torsional stiffness 


mass moment of inertia of fixed weight, about CG 

"s 

mass moment of inertia of store about pivot point 

"a 

mass moment of inertia of wing about E.A., per 
unit length 

k 

reduced frequency bw/U 

1 

wing semi- span 


spanwise location of store , based on 1 

''p 

chordwise location of fixed weight c.g., measured 
from EA, based on b, and positive aft 

's 

chordwise location of store pitch pivot, measured 
from wing mid-chord, positive aft, based on b 


vertical distance of fixed weight, measured from 
wing mean chord line, based on b, positive dowxi 


vertical distance of store pitch pivot point, 
measured from mean chord line, based on b, positive 
down 



3 


m 

IHj 

m. 


’a 


a 


mass of wing per unit length 
mass of fixed weight 
mass of store 

wing radius of gyration based on wing half chord 
and referred to wing EA 

store radius of gyration, based on wing half chord, 
referred to store pitch axis location 

radius of gyration of fixed weight, based on half 
chord, referred to wing E.A. (-mpb^/I^^) 


forward velocity 


wing center of gravity location, based on wing 
half chord, positive aft 


X 


s 


store center of gravity location, based on wing 
half chord, positive aft 


X 

K 

P 


Wj^(N) 

%s 

“ts 

( 1 ) 


spanwise coordinate, based on I 

2 

virtual mass of wing Trpb /m 
air density 

fundamental uncoupled banding frequency of wing 
with rigid mass 

bare wing uncoupling bending frequency in mode 

'th 

bare wing uncoupled torsion frequency in M mode 
uncoupled store vertical frequency 
uncoupled store fore and aft frequency 
uncoupled store pitch frequency 
flutter frequency 



Concept Description 


The concept of installing a flexibly mounted store on an 
aircraft to increase the flutter speed is based mainly on the 
knowledge that flutter is a function of the torsional frequency 
In general, the higher the torsional frequency, the higher the 
flutter speed. When a large store is rigidly fastened to a 
wing, the torsional frequency is reduced which results in a 
reduction in flutter speed. If a system of store suspension 
could be devised such that the store inertia effect was 
essentially removed from the system, it was presumed that the 
flutter speed should show a marked increase since the torsmnal 
frequency would not suffer a large reduction. A system to 
accomplish this is mechanically very easy, for it is only 
necessary to hinge the store on a pivot with low flexibility. 
Fig. 1 illustrates the mechanical arrangement of such a system. 

From practical considerations, it would not be advisable 
to maintain a store in a completely free condition, so some 
spring stiffness is required to maintain the store in its 
proper relationship to the wing. Also some system for maintain 
ing the alignment of the store with the wing may be necessary. 
Such a system should have a very low frequency in order to 
maintain the low pitch frequency of the store. The present 
report analytically investigates this concept and presents some 


of the results . 
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Analysis 

The analysis was specialized to the case of a structurally 
uniform wing, having a rectangular aerodynamic planform of 
relatively high aspect ratio. Three degrees of freedom in wing 
torsion, and four degrees of freedom in wing bending were 
utilized in a modal type analysis. A single, spring mounted 
weight was added to the system which could be located at any 
spanwise or chordwise position. The spring-mass was restricted 
to three degrees-of-freedom, namely, pitch, vertical translation, 
and fore-and aft motion. Two-dimensional incompressible, 
unsteady aerodyneimic theory was used, ref. 1. Details of the 
analysis are given in the Appendix. The next section contains 
a discussion of the results. 

Results and Discussion 

The analysis technique presented was applied to a rec- 
tangular wing whose properties were held constant throughout 
the study. The wing characteristics are given in Table 1. 

A small store support bracket was assvimed to be rigidly attached 
to the wing and its properties are given in Table 2. The 
store properties are given in Table 3 . 

Of primary importance is the store pitch stiffness and 
this property was first considered, the results of which are 
given in Figure 2. The variation of center of gravity position 
is shown in Figure 3, and finally the effect of moving the 
store spanwise is given in Figure 4. 
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Effect of Store Pitch Frequency on Flutter 

On Figure 2, the ratio of the flutter dynamic pressure 
for the wing with a flexible store to the flutter dynamic 
pressure of the wing with no store attached is plotted against 
the ratio of the store uncoupled pitch frequency to the wing 
bending frequency with the rigid store attached. Starting 
at the left of the figure, tOQ/a)j^=0, it is seen that for a 
completely free store in pitch, the flutter speed is well 
above the bare wing flutter speed. As the store pitch spring 
stiffness is increased, the flutter speed rapidly increases 
to until (jOq/cOj^-I; for largei^values of Wg/tOj^, a rapid decrease 
is found which has a minimum at 1*5, where the flutter 

speed begins to slowly increase and approaches an asymptotic 
value for tOo/w, - which represents a rigid pylon. There- 
fore, for this configuration, if the ratio is maintained 

below unity, substantial benefits arise from the use of a low 
store pitch frequency. It should be pointed out that it was 
found that the pylon stiffness in the vertical direction must 
be stiff, since some calculations made with a reduced vertical 
translation stiffness reduced the flutter speed well below the 
base wing flutter speed. 

Effect of Store Center of Gravity Location 

In Figure 3 is shown the results of calculations for 
a rigid and a reduced stiffness case shown in Figure 2, 
except that the center of gravity of the store was moved 
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fore and aft. Examining the rigid pylon case first, it is 

seen that as the C. G. is moved forward from an aft position, 

the flutter speed is essentially constant in the range of 

■ of from .3 to .5, which is well below the flutter 
base wing 

speed of the bare wing. As the center of gravity approaches 

a more forward position, the flutter speed increases rapidly 

when X == -. 1 . 
s 

For the case of reduced stiffness in the store pitch mode 

(Wg/Wj^ = 0.54 when x^ = 0) the flutter speed trend is virtually 

independent of C. G. location, ranging from wing 

for X = .4 to q/q. . =1.7 for x = -.4. 

s ^ ^bare wing s 

Effect of Spanwise Store Location 

Calculations were made to determine the effect of the store 
spanwise location for both the rigid and flexible pylon. These 
results are given in Figure 4. For the rigid pylon, the flutter 
q drops as the store is moved from the base to the wing midspan, 
then increases as the wing tip is approached and for the store 
located at the wing tip, the bare wing flutter speed is approached. 
For the reduced stiffness pitch case, (“ 3 / 03 ^^ = 0.54) the flutter 
speed gradually increases as the store is moved towards the wing 
tip where the flutter q is about 2.4 times the bare wing flutter q. 

Concluding Remarks 

A calculation procedure has been presented for a uniform 
cantilever wing on which is attached a flexibly mounted store. 
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Results of calculations based on the procedure indicated the 
following results: 

1) The flutter speed is well above the bare wing 

flutter speed for frequency ratio less than 

unity. 

2) Over the range of pylon stiffnesses, representative 
of production aircraft, the flutter speed was con- 
siderably less than that of the bare wing. 

3) Large variation in some of the more important store 
parameters showed the flutter speed to be virtually 
independent of the variations and the speed remained 
well above the bare wing flutter speed, if the ratio 
a)„/co, is maintained 

9 n 


below unity. 
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APPENDIX A 
Equations of Motion 

The equations of motion for flutter were derived on the 
basis of Lagrange's equations. In matrix form, the equations 
appear as A = XB where the coefficients A & B are given in detail 
subsequently, and X=(j^ ) (1+ig) . The equations were specialized 

for the case of a uniform cantilever wing carrying a fixed and 
sprung mass at an arbitrary spanwise and chordwise location 
(figure 1) . The equations were linearized and the unsteady 
aerodynamics were limited to the incompressible, two dimensional 
case. For the structural part, ten degrees of freedom were 
utilized, namely, four uncoupled uniform beam bending modes, 
three uncoupled uniform beam torsion modes, and three degrees 
of freedom for the sprung mass comprising pitch, heave, and 
fore-and aft degrees of freedom. No aerodyneunic forces were 
applied to the store. 

Since the wing mass, inertial and unsteady aerodynamic 
forces are constant along the span, and theoretical uncoupled 
beam modes were assumed, it was convenient to calculate the 
generalized mass, inertial and aerodynamic terms in closed form. 
The expressions used for the uniform cantilever modes were 
slightly different from those found in textbooks on vibration. 

It has been found that the forms given below provide more accurate 
results for the higher modes. 
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For odd modes, the bending modal function is 


^ cosh(B /2)(2x-l) sin(6^/2)(2x-l) 

= 2 ^ cosh(B /2) sin(B„/2) ^ 


For even modes , the modal function is 


^ sinh (By2X2x-l) cos (By.2)(2x-1) 

" 2^ sinh B /2 cos 6„/2 ^ 


n 


n' 


where B^ must satisfy the following relation 


cosh B cos B +1=0 n=l,2,3,4... 

n n 


and where x is the spanwise coordinate, normalized by the span 
length , i . 

For the torsion mode, the modal function was selected to be 
fg(M) = (-D^'^^sin ir(m-i) x m=l,2,3... 


Several integrals of these functions are required as 
follows ; 

For N odd 


FW(N) = dx 

0 

(sinh B_,/2) cosh B„/2+B^/2 

ii II n 

4 B^ Cosh‘S B^/2 




6 /2~ sin B_ 
ii n. 

4B sin“ B_/2 
n n 
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For N even 

1 2 

FW(N) = /q f^(N) dx 

sinh (Bj^/2)cosh 

43^ sinh^ B^/2 

■=■ sin B„ 

+ 2 n 

43 cos^ 3„/2 
n n 

For N odd 

1 

F(N,M) = / f (N)f (M) dx 
0 

3 cotan(3 /2) - (m-1/2) 

_ n n 

-23 ^ + 77 ( 111 - 1 / 2 )^ 

n 

B tank(3 /2 ) +tt (m-2/2) 

+ n n 

23 ^ + 277 (171-1/2)^ 
n 


For N even 

1 

F(N,M) = / f„(N) f (M) dx 

0 

3 cotanh ( 3 _/ 2 ) - (-1) ^^^^^ni-1/2) tt 
_ n n 

2(3jj^ + 77(m-l/2)^^ 

1-3 tan3„/2+ (-1 ) "'‘^^77 (m- 1/2 ) 

^ n n 

2(-3j^^ + 77 (m- 1/2)^) 
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The order of the inodes in the matrix is as follows - the 
first four rows are the bending modes / the 5—7 rows are the 
torsion modes, and the 8-10 rows are store vertical translation, 
store pitch, and store fore and aft respectively. 

The definition of the elements of the matrices A and B 
follows where the symbols are defined after the definition: 
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MATRIX DIAGONAL ELEMENTS 
N=M N= 1 , 2 , 3 , 4 

A(N,M) = -FW(N) (Hl)-(mgAJ (FWL(N) + (-mpAn^) (FWL (N) ) ^ 
+ i (2KF(FW(N) )/k) 

N=M N= 5,6,7 

A{N,M) = -FA(N-4) (H12) - Rg (FAL (N-N) ) 

+ {FI) (FAL(N-4) + i(FA(N-4) (H13) ) 


A(8,8) = -1 


A(9,9) = FAL(l) 

A(10,10) = -1 

OFF DIAGONAL TERMS 

N= 1,2,3, 4 M=l,2,3,4 M?^N 

A (N,M) = - ( (mp4mg) /H^) (FWL (N) ) (FWL (M) ) 

N=5,6,7 M=5,6,7 M?^ 

A(N,M) = (-(m^/il„)r„ + FI) (FAL(N-4) (FAL(m-4)) 

S lU ij 
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N=5,6,7 M=l,2,3,4 

A(N,M) = (F(M,N-4) ) (HIO) 

- (H7) (FAL(N-4) ) (FWL (M) ) /r^ 

+rj FWL(M) - FAL(M-4) 

r r 

+i(- (F(M,N-y) ) (Hll) 


N=1 ,2,3,4 M=5 , 6 , 7 

A(N,M)=F(N,M-4) (H5)-(H7) (FAL(l)) (FWL(M-4) 

- (M„/£„(FAL (N) ) (FWL(M-4)) 
r r 

+ I (F(N,M-4) ) (H6) 


N=l,2,3,4 M=8 

A(N,M) = - {M^/£„) (FWL(N) ) 

S A’l 


N=l,2,3,4 M=9 

A(N,M) = -(Mg/S-j^)Xg(FWL(N)) 


A(l,10) = A(2,10) = A(3,10) = A(4,10) = 0 


N=5,6,7 M=8 

A(N,8) = -(H7) (FAL(N-4)/r^ 


N=5,6,7 M=9 

A(N,9) = (H16) (FAL(N-4) ) 
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N=5,6,7 M=10 

A(N,10) = (FAL(N-4))/r^ 


N=8 M=l,2,3,4 

A(8,M) = FWL(M) 


N=8 M=5,6,7 

A(8,M) = - (£^+x^) (FAL(M-4) ) 
s s 

A/8,9 = -Xg 


N=9 M=l,2,3,4 

A(9,M) = (H14) (FWL(N)) 


N=9 M=5,6,7 

A(9,M) = (H15) (FAL(M-4)) 


A(9,8) = H14 


A{9,9) = -1 


A(10,l) = Ado, 2) = A(10,3) = Ado, 4) = 0 


N=10 M=5,6,7 

A(10,M) = £j^(FAL(M-D) 


A(10,10) = -1 


N=1,2,3,4 


2 


B(N,M) = -( 




(N) 


u 


a. 


(FW(N) ) (1 + ig (N) ) 


N=5 

B(5,5) = -FA(1) (1 + ig(5) ) 


B(N,N) = - ) (FA(N-4))( 1 + ig(N)) 



B(8,8) = - (^) ( 1 + ig(8)) 
“l 


2 

B(9,9) = -(— ) (1 + ig (9)) 

^1 


2 

“t 

B(10,10) = - (y— ) ( 1 + ig(lO)) 
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HI = 

H2 = 
H3 = 

H5 = 

H6 = 
H7 = 
H8 = 

H9 = 

HIO 

Hll 

H12 

H13 

H14 

H15 

H16 



FI = 


2k 

1 + K + 


2k „ 

ir ^ 


m /Jl 
s m 


-X + Ka + - 2KG(i - a)/k 

QL 1 ^ ^ 


K/k + 2<G/k^ + 2KF(i - a)/k 


(Z + X ) m /Z 
s s s m 

/ 2 


2y 2 

X /r 
s' s 


-x^/r^ +Ka/r^ + 2K(a+i) G/(kr^) 


2K(a+|) F/(kr^) 

1 + K{^a^)/rJ + 2 {9+|)F/(kV) 

2K(a+i) (a“)G/(kr^) 

(k(| - a)/k - 2k (a + |)G/k^ - 2KF(a + |) (a 
-"s/^s 


|)/k)/r2 






2 2 2 2 
Z7 + X + Z + Z X + r 
1 s s s s s 


-((&F + )iiipbV(V) + Ip/(Ij^£)) 
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TABLE 1. - WING CHARACTERISTICS 

X = .15754 
a 

K = .06204 

10 2 

El = 2.756 X 10 N/m 
GJ = 2.407 X 10^ N/m^ 
m =59.13 Kg/m 
I = .00282 Kg-m^ 

= .1549 
a = - . 3 

TABLE 2. - FIXED WEIGHT CHARACTERISTICS 

.5922 Kg 
.1473 

0 

2 

.00727 Kg-m 

TABLE 3. - STORE CHARACTERISTICS 

x ^ = -.01 

S 

r? = .4162 
s 

m = .0382 Kg 
s 

I = .2877 Kg-m^ 
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Figure 2.- 


Effect of store pitch frequency on wing-store flutter. 
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Effect of spanwise stores location. 


Figure 4 . 





